Stretching of a polymer below the point 
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The unfolding of a polymer below the point when pulled by an external force is studied both 
in d = 2 on the lattice and in d = 3 off lattice. A ground state analysis of finite length chains shows 
that the globule unfolds via multiple steps, corresponding to transitions between different minima, 
in both cases. In the infinite length limit, these intermediate minima have a qualitative effect only in 
d — 2. The phase diagram in d — 2 is determined using transfer matrix techniques. Energy-entropy 
and renormalization group arguments are given which predict a qualitatively correct phase diagram 
and a change of the order of the transition from d — 2 to d — 3. 



The recent refinements in experimental techniques em- 
ploying optical tweezers Q , atomic force microscopes [[| , 
soft microneedles jjj, make potentially possible for re- 
searchers to monitor the behaviour under tension and 
stress of various biopolymers and then to elucidate the 
mechanism of some force-driven phase transitions occur- 
ring at the single molecule level, such as the unfolding of 
the giant titine protein B, the stretching of single col- 
lapsed DN A molecules , the unzipping of RN A and 
DNA Theoretically, on the other hand, quite a few 
statistical mechanics models have been subsequently pro- 
posed in order to explain these experimental results and 
to identify the physical mechanisms behind these phase 
transitions [THlq]. In this work we focus on two theo- 
retical models aimed at understanding the unfolding be- 
haviour of polymers in a poor [|l6| solvent (globules), i.e. 
below the theta temperature, Tq. 

This is motivated by the vast number of experimen- 
tal results in the literature. In particular, various force 
vs elongation (/ vs x) curves have been recorded in 
experiments studying these phenomena ||||. In most 
cases a curve consisting of three distinct regimes, and in 
articular displaying a plateau for intermediate stretch 
, have been observed; whereas in a few examples, for 
shorter polymers, a stick-release pattern with histeresis 
has been found. The first observation is in good agree- 
ment with the mean field theory proposed in Ref. , and 
the plateau strongly suggests the presence of a first or- 
der phase transition (see also the recent measurements in 
||). Our results suggest that there might be more than 
one possible shape for the / vs. x curves according to d, 
the spatial dimension and to polymer length, so that a 
different behaviour occurs when mean field theory @ is 
qualitatively incorrect. These models have a remarkable 
interest even on a purely theoretical ground. First, the 
numerical study recently performed in has suggested 
the possibility that the transition is second order in d = 2 
and first order in d = 3. This has been confirmed to some 
extent in a study of the model on a hierarchical lattice 
with fractal dimension two . The d = 2 case is impor- 
tant as it is below the upper critical dimension for theta 
collapse and mean field predictions may well be incorrect. 
A thourough analysis and a clear physical mechanism un- 



derlying the difference of the nature of the transition as 
d changes are then needed. Second, the mean field anal- 
ysis of Ref. ]l0| has suggested there could be a re-entrant 
region in the phase diagram for low temperature similar 
to what happens for DNA unzipping Jl5[] . However the 
exact results in prove mean field is not valid in d = 2. 

Then, our aim here is to describe theoretically the un- 
folding transition of globules not relying on the mean 
field approximation. We first characterize the evolution 
of the ground states of a finite polymer as the pulling 
force increases. This will help to understand the thermo- 
dynamics. In particular, we compute the phase diagram 
in the temperature-force (T, /)-plane in d = 2 on the lat- 
tice, where we can use the transfer matrix method (Fig. 
2) together with exact enumerations (Fig. 1). 

We begin by considering a self- avoiding walk (SAW) on 
the square lattice with fixed origin. The model partition 
function (generating function) in the canonical ensemble 
in which T and /, the stretching force, are fixed is: 
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where TV is the number of monomers (including the ori- 
gin) of the SAW, H(C) (referred to a configuration C, 
i.e. to a SAW) is the energy of a SAW, b is the number 
of pairs of neighboring occupied sites not adiacent along 
the chain, (3 = ^ is the inverse temperature, w n(N, b) is 
the number of configurations of a SAW with fixed origin 
o and end-to-end distance R, of lenght N and b contacts, 
R x is the projection along the force direction (x axis) of 
R. We take both the Boltzmann constant and e equal to 
one. 

When T is low, one may look for the ground states 
among the rectangles of sides L x and L y that are com- 
pletely covered by the SAW (in other words such that 
L x L y = N, where we neglect the small effects arising 
when this rectangle cannot be constructed with both 
L XjV integers). The energy of this rectangular hamilto- 
nian walk with a non zero / is —H(N = L x L yi L x ) = 
N - L x - £- + 1 + f(L x - 1). The minimum of H(N, L x ) 
for given N with respect to L x yields the most stable 
configuration for various values of /. The minimum oc- 
curs for an /-dependent value of l? x = L®(N, /), namely 
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L X (N, f) = yy^-j. For any / < 1 one has a compact 

configuration. However, when the critical value / = 1 
(for T = 0) is reached all integer values of L x from N 
(stretched coil) to L x ~ TV 1 / 2 (compact globule) become 
degenerate for large N . Note that this does not hold in 
d = 3 where it is well known that there is a Rayleigh 
instability in the thermodynamic limit P,|7|,pT| . This can 
be seen by comparing the globule energy - which in d = 3 
is 2N in the large TV limit - with the energy of a paral- 
lelepiped, with elongation along / equal to L x and with 
edges L y = L z in the perpendicular plane. The force 
above which (for T = 0) the parallepiped is a better 
ground state than the three-dimensional globule is 2L y 
(L y <C A/ -1 / 3 ), and there is no longer any degeneracy at 
the critical force f c = 2 (at T = 0). 

In Fig. |l|a we sketch the situation in d = 2. The 
minima hierarchy, shown in the shaded area in the top 
panel, affects the low T region of the < x > (average 
elongation) vs. / curves for finite length (bottom panel). 
However, only one transition survives in the large N limit 
and represents a true phase transition (as represented in 
Fig. [l]a by the shaded wedge ending in just one point 
in the TV -1 = axis). Similarly, when T is raised the 
multi-step character of the < x > vs. / curves is lost due 
to fluctuations which blur the ground state dominance in 
the partition. In Fig. [l]b we show the analogous picture 
for a 3d— model discussed below. Short chains in d = 2 
and 3 behave similarly whereas in the infinite length limit 
the d = 3-case shows an abrupt unfolding transition. 

We now turn to the thermodynamic behaviour of the 
SAW model on a square lattice. We use the transfer ma- 
trix (TM) technique, following Refs. 

Let us introduce briefly the principal features of the 
TM approach: the partition function of a polymer of 
N sites is given by Eq. (|l|), with / = 0. In the ther- 
modynamic limit (N — > oo ) we expect that Zn(T) ~ 
\p(T)] N , then the free energy per monomer T is simply, 
T = — T log /i(T). It is more convenient [[l9| to introduce 
the following generating function (z is the step fugacity) 

SoR = X> J V 6 «; oR (JV ) &) (2) 
N,b 

It is known that for z < z c (T) = l//i(T), the in- 
verse SAW connectivity, # r ~ exp(—R x /£(z, T)), where 
£(z,T) is the correlation length and R x is the projection 
of R along x. We study the stretching of an interacting 
SAW in a strip of finite size L along y and infinite length 
along x. It is possible to define [|19[ an L-dependent 
correlation length ^l(z^T) via the formula £l(z,T) = 
— \ oz Xl(z t) •> wnere ^l(z,T) is the largest eigenvalue of 
the transfer matrix, that equals 1 at z = z£(T). We 
apply the phenomenological renormalization, to find suc- 
cessive estimates for z c (T) = lirn^oo z£(T). Including 
the force via Eq. (|l|), the equation for the critical force 
f c (T) is then ideally found via: 

fc — —T lim lim logAz,. (3) 

z^z c (T)- L^oo 



The order of the limits in Eq. (||) and a correct choice of 
the boundary conditions (see below) are crucial |2^j . 

In Fig. H the phase diagram for the stretched inter- 
acting SAW is shown. With the TM, a right choice 
of the boundary conditions is needed |]19|-pi|. We have 
used both periodic (PBC) and free boundary conditions 
(FBC). PBC have been employed to get the best es- 
timate of z c (T) through phenomenological renormal- 
ization. This value is then used with FBC, to find 
the correct L-dependent critical eigenvalue \l(z c {T),T). 
Finally, adopting the extrapolation algorithm of [[23[, 
lmiL^oo \l(z c (T),T) is obtained which, through Eq. 
allows to get the phase diagram. As expected, in the case 
of FBC, there are oscillations in data going from odd L to 
even L. As usual in this context, a separated analysis of 
even and odd L data was necessary for obtaining a bet- 
ter convergence (see Fig. ^|). One point on the transition 
line obtained previously in jnj, is recovered here. 

One can get an approximate description of the tran- 
sition if one requires that the globule and coil phases 
coexist. The globule free energy is easily estimated 
in terms of hamiltonian walks fig] . On a square lat- 
tice the energy is simply given by minus the length of 
the polymer whereas the entropy is given in terms of 
the number of hamiltonian walks which grows exponen- 
tially with N fl6|| . Thus the globule free energy per 
monomer is T g = —1 — Tlog(4/e) where we have used 
the accurate mean field estimate of the entropy as given 
in [ p4| . The coil free energy T c is approximated as 
that of an unconstrained random walk in presence of a 
pulling force and contacts are neclected. We thus get 
T c = — log(2(l + cosh(/3/))). At coexistence one finds 
f c (T) = Tcosli (2 exp(l/T— 1) — 1) (the continuos curve 
in Fig. §. 

We note that / c (0) = 1 is the exact result and at low 
T the phase diagram displays a reentrant region. As 
T — > T^i, / approaches rather smoothly, consistently 

with the prediction f ~ (T e - T)^ = ^ @. Within the 
TM approach one can also infer the order of the transi- 
tion. To do this we observe that, if for z — > z c (T)~: 

X(z c (T), T) - X(z, T) ~ (z c (T) - z) A , (4) 

then A < 1 (A = 1) means a second (first) order transi- 
tion, as < x >~ (/ — /c)^ -1 for / ~ f c . Our data at not 
too low T are compatible with a second order transition 
(inset of Fig. ||). 

Inspired by exact renormalization group (RG) on the 
Sierpinski lattices O], and on approximate RG in d— 
dimensional lattices|2q1 , we propose the following sim- 
plified real space RG which rationalizes our results. The 
RG recursions relations can be written for the generating 
function representing polymers traversing a hypercube of 
linear size 1 once, A, and 2 d_1 -times, B (Fig. 3). The 
terms A and B represent parts of the chain which are 
in the coil and globular state respectively. The recursion 
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relations can be calculated as in Ref. or by enumer- 
ating the SAWs on 2x2 or 2x2x2 cells as in j25|. To 
leading order in A at ^ T <C T e 

A' = A 2 + a(d)A 2 B 2 , B' = j(d)B 2 \ (5) 

where a(d) and 7(d) are d— dependent constants. There 
are three fixed points in the flux in Eq. ||: A — 0, 

B = B* = 7(d) 2d ~ 1 corresponds to the globular phase, 
(A, B) = (1, 0) to the coil phase, while the last non-trivial 
fixed point A = A* = i +a ^^*) 2 ? ^ ~ ^* cnarac terizes 
the unfolding transition. The value of the term a(d) af- 
fects the behaviour of the RSRG flux near the fixed point 
(A*, 5*). One can see {2£§ that a(d) ^ if d = 2 and is 
in d > 2. When a(d) ^ (i.e. in d = 2), the RG flux is 
smooth near (A*, 5*) and the critical fugacity near / = 
f c (T) behaves as z c (f c (T),T) - z c (f,T) ~ (/ - / C (T)) 2 
signalling a second order transition with A = 1/2 in Eqfl. 
On the other hand when a(d) = (i.e. d > 2), the tran- 
sition is first order and two-state like. The presence of 
the mixed term in Eq. (||) is crucial and enhances the 
entropy of the coil phase since it contributes to A! . Con- 
sequently, the d = 2 two-state approximation in Fig. || 
gives a transition line which is higher than the numerical 
result for ^ T <C Tq. The entropy gain in the stretched 
coil, as T — > 0, is hampered as it costs a finite surface 
energy (dominant as T — > 0) to change locally an elon- 
gated globular region into a coil and vice-versa. This is 
why the solid curve in Fig. 2 matches our numerics as 
T 0. 

Let us now discuss the d = 3 case. The model we used 
is the freely jointed chain (FJC) (see e.g. |u|) in the con- 
tinuum (off lattice). The FJC is subject to a compact- 
ing pairwise attractive potential between non-consecutive 
beads and to a stretching force / at the extrema. The 
pairwise potential is chosen to be an asymmetric square 
well with a hard core radius, 2i?^ c , which acts between 
non-consecutive beads along the chain, and an attraction 
range Ri > 2Rh c , i.e. the distance up to which the parti- 
cles interact. We have checked that the results reported 
in the following do not appreciably depend on the two- 
body potential details. The parameters we have used to 
generate the configurations shown in Fig. |l|b (top panel) 
are Rh c — 0.6, R\ = 1.6 where the unit length is the 
distance between sucessive beads along the chain. 

The ground states of short chains (up to N = 30 bead 
long) has been determined by performing simulated an- 
nealing employing Monte-Carlo dynamical simulations. 
The FJC is evolved dynamically by means of three sets 
of moves: the pivot, reptation and crankshaft moves [p7| . 
We lower T during the simulation according to a standard 
annealing schedule. We found evidence also in this d = 3 
case that the unfolding of a finite length homopolymer 
proceeds in a multistep fashion. The collapsed globule 
first orients itself along the pulling force as soon as there 
is a nonzero /. At larger / the globule is slightly elon- 
gated (much less than in the d = 2 case) and after this 
a helix forms followed by a zig-zag curve and finally by 



a stretched coil. This succession of minima, shown in 
fig. 0b is intriguing for a two-fold reason: firstly because 
it suggests that helices, one of the well known building 
block of proteins, come out rather naturally as one of 
a few stable minima of a homopolymer in a poor sol- 
vent subject to a finite stretching force; secondly because 
the unfolding transition of a finite polymer in this model 
appears to be markedly different from a globule-to-coil 
(two-state) transition. The mean field picture of an all or 
none transition is recovered for infinite length. The situ- 
ation is depicted in Fig. |l|b. Since in d = 3 the transition 
is first order, the mean field treatment is correct in the 
thermodynamic limit. The fact that helices become bet- 
ter ground states than compact globules could be easily 
verifiable in experiments done in the fixed force ensemble. 
These can now be confidently performed, with not more 
substantial difficulties than the more conventional ones, 
performed in the fixed stretch ensemble. Indeed the stick- 
release pattern reported in Ref. ^ shows several different 
peaks in the / vs. < x > curve, and is compatible with 
a 'multi-step transition' in which the globule undergoes 
more than one conformational changes (signalled by the 
peaks) during the unfolding. Intriguingly, this pattern is 
reported for smaller polymers, a fact which would be in 
agreements with our findings; whereas for large N this 
effect is much less important and the curves have a single 
(dominant) plateau. 

Helices appear as ground states for a potential con- 
sisting of a force term and a hydrophobic contact poten- 
tial because they are both elongated and offer a good 
shielding from the outside solvent to the monomers ly- 
ing in its interior. Note also that helices arise as opti- 
mal states of tubes of non-zero thickness subject to com- 
paction [ p8[ . If a worm-like chain is used instead of the 
FJC, the picture should not change, because the elon- 
gated states such as the helices and the zig-zag are stiffer 
than the collapsed globule (sharp corners between succes- 
sive monomers are necessary in this state) and so should 
be even more favoured as / is increased. 

In conclusion, we have discussed the unfolding transi- 
tion of a homopolymer under the action of an external 
pulling force in d = 2 on the lattice and in d = 3 off lat- 
tice. A ground state analysis shows that for finite length 
polymers, the unfolding is not abrupt, rather it occurs 
via a multi step sequence of states. These are more elon- 
gated than the globule and make more contacts than a 
coil. In d = 3 helices arise naturally as ground states at 
intermediate forces. For infinite polymers, on the other 
hand, the situation is different: in d = 3 the intermedi- 
ate ground states disappear due to the Rayleigh instabil- 
ity, and the transition is effectively two-state, whereas in 
d = 2 they survive in the thermodynamic limit. Indeed, 
from the transfer matrix results, it is apparent that the 
mean field hypothesis is incorrect in d = 2 even at rather 
low T: the transition is second order as also found in 
Ref. (IT) although the prediction of the re-entrant region 
agrees with the TM results. Furthermore a renormal- 
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ization group based argument is in agreement with this behaviour. Details will be given in a future work. 
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FIG. 1. a) d — 2 Schematic diagram of minima for poly- 
mers of different sizes (top) and low T < x > vs. / curve for 
N = 25 (bottom) found with exact enumerations in d = 2. 
b) d=3: Same as in a), except that the < x > vs. / curve 
(bottom) is for N = 30 and is found by simulated annealing. 
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FIG. 2. Phase diagram for the stretching of a SAW on a 
2d- lattice, obtained with the TM technique. Inset: plot of 
log A vs. log z for T = 0.7. 
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FIG. 3. Representation of the generating functions A and 
B in the text (for d = 3) at a generic level of iteration in 
the RG. Thick lines represent portions of the SAW that enter 
through a vertex and go out from another one. 
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